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Abstract 

The nucleus of 208 Pb — a system that is 18 order of magnitudes smaller and 55 orders of 
magnitude lighter than a neutron star — may be used as a miniature surrogate to establish im- 
portant correlations between its neutron skin and several neutron-star properties. Indeed, a nearly 
model-independent correlation develops between the neutron skin of 208 Pb and the liquid-to-solid 
transition density in a neutron star. Further, we illustrate how a measurement of the neutron 
skin in 208 Pb may be used to place important constraints on the cooling mechanism operating in 
neutron stars and may help elucidate the existence of quarks stars. 
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I. INTRODUCTION 



It is an extrapolation of 18 orders of magnitude from the neutron radius of a heavy nucleus 
- such as 208 Pb with a neutron radius of _R n ~5.7 fm — to the approximately 10 km radius 
of a neutron star. Yet both radii depend on our incomplete knowledge of the equation of 
state of neutron-rich matter. That strong correlations arise among objects of such disparate 
sizes is not difficult to understand. Heavy nuclei develop a neutron-rich skin as a result of 
its large neutron excess (e.g., N/Z = 1.54 in 208 Pb) and because the large Coulomb barrier 
reduces the proton density at the surface of the nucleus. Thus the thickness of the neutron 
skin depends on the pressure that pushes neutrons out against surface tension. As a result, 
the greater the pressure, the thicker the neutron skin [l|. Yet it is this same pressure that 
supports a neutron star against gravitational collapse 0, El • Thus models with thicker 
neutron skins often produce neutron stars with larger radii |4j. 

The above discussion suggests that an accurate and model-independent measurement of 
the neutron skin of even a single heavy nucleus may have important implications for neutron- 
star properties. Attempts at mapping the neutron distribution have traditionally relied on 
strongly-interacting probes. While highly mature and successful, it is unlikely that the 
hadronic program will ever attain the precision status that the electroweak program enjoys. 
This is due to the large and controversial uncertainties in the reaction mechanism p, 
The mismatch in our knowledge of the proton radius in 208 Pb relative to that of the neutron 
radius provides a striking example of the current situation: while the charge radius of 208 Pb 
is known to better than 0.001 fm 0], realistic estimates place the uncertainty in the neutron 
radius at about 0.2 fm 0]. 

The enormously successful parity-violating program at the Jefferson Laboratory j^, 
provides an attractive electroweak alternative to the hadronic program. Indeed, the Parity 
Radius Experiment (PREX) at the Jefferson Laboratory aims to measure the neutron ra- 
dius of 208 Pb accurately (to within 0.05 fm) and model independently via parity-violating 
electron scattering Parity violation at low momentum transfers is particularly sensi- 
tive to the neutron density because the Z° boson couples primarily to neutrons. Moreover, 
the parity-violating asymmetry, while small, can be interpreted with as much confidence 
as conventional electromagnetic scattering experiments. PREX will provide a unique ob- 
servational constraint on the thickness of the neutron skin of a heavy nucleus. We note 
that since first proposed in 1999, many of the technical difficulties intrinsic to such a chal- 
lenging experiment have been met. For example, during the recent activity at the Hall A 
Proton Parity Ex per iment (HAPPEX), significant progress was made in controlling helicity 
correlated errors [llj. Other technical problems are currently being solved — such as the 
designed of a new septum magnet — and a specific timeline has been provided to solve all 
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remaining problems within the next two years 

Our aim in this contribution is to report on some of our recent results that examine the 
correlation between the neutron skin of 208 Pb and various neutron-star properties 0, 2, 13 ■ 
In particular, we examine the consequences of a "softer" equation of state that is based on a 
new accurately calibrated relativistic parameter set that has been constrained by both the 
ground state properties of finite nuclei and their linear response. Further, results obtained 
with this new parameter set — dubbed "FSUGold " |l4l| — will be compared against the NL3 
parameter set of Lalazissis, Konig, and Ring [3, Il6l| that, while highly successful, predicts 
a significantly stiffer equation of state. 
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Model 


m s gl gl g 2 p k X ( A v 


NL3 


508.1940 104.3871 165.5854 79.6000 3.8599 -0.0159 0.0000 0.0000 


FSUGold 


491.5000 112.1996 204.5469 138.4701 1.4203 +0.0238 0.0600 0.0300 



TABLE I: Model parameters used in the calculations. The parameter k and the inverse scalar 
range m s are given in MeV. The nucleon, omega, and rho masses are kept fixed at M = 939 MeV, 
m u = 782.5 MeV, and m p = 763 MeV, respectively. 



II. FORMALISM 

The starting point for the calculation of the properties of finite nuclei and neutron stars 
is an effective field-theory model based on the following Lagrangian density: 



Ant = V> 
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-^f-^Y+^glv^) +K(gl\.w)(glv p y») . (1) 

The Lagrangian density includes an isodoublet nucleon field (ip) interacting via the exchange 
of two isoscalar mesons — a scalar (0) and a vector (V*) — one isovector meson and the 
photon (v4 M ) [ttI. Ei|. In addition to meson-nucleon interactions, the Lagrangian density is 
supplemented by four nonlinear meson interactions, with coupling constants denoted by k, A, 
(, and A v . The first three of these terms are responsible for a softening of the equation of state 
of symmetric nuclear matter at both normal and high densities [HI. In particular, the cubic 
(k) and quartic (A) scalar self-energy terms are needed to reduce the compression modulus of 
symmetric nuclear matter, in accordance to measurements of the giant monopole resonance 
in medium to heavy nuclei j^]. In turn, tu-meson self-interactions (£) are instrumental for 
the softening of the equation of state at high density — thereby affecting primarily the 
limiting masses of neutron stars [3]. Finally, the last of the coupling constants (A v ) induces 
isoscalar-isovector mixingand has been added to tune the poorly-known density dependence 
of the symmetry energy 0, I3| • As a result of the strong correlation between the neutron 



radius of heavy nuclei and the pressure of neutron-rich matter [jj, |21[ , the neutron skin of a 
heavy nucleus is highly sensitive to changes in A v . 



III. RESULTS 

In Table H] we list the various mass-parameters and coupling constants for both of the 
models (NL3 and FSUGold) employed in the text. Further, in TablelTDa comparison is made 
between the very successful NL3 parametrization|l5l Il6j|. FSUGold, and (when available) 
experimental data. Finally, in Fig. ^we display the charge and neutron density of 208 Pb as 
a function of the isoscalar-isovector parameter A v . Note that while A v modifies the neutron 
radius of a heavy nucleus, it does so without affecting those ground-state properties that 
are well constrained by experiment, such as binding energies and charge radii. 

While the agreement (at the 1% level or better) between NL3 and experiment is sat- 
isfactory — and this agreement extends all over the periodic table - incorporating 
additional terms into the density functional is demanded, not by ground state data, but 
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Nucleus 


Observable 


Experiment 


NL3 


FSUGold 


40 Ca 


B/A (MeV) 


8.55 


8.54 


8.54 




R ch (fm) 


3.45 


3.46 


3.42 




Rn-Rp (fm) 


— 


-0.05 


-0.05 


48 Ca 


B/A (MeV) 


8.67 


8.64 


8.58 




R ch (fm) 


3.45 


3.46 


3.45 




R n -Rp (fm) 


— 


0.23 


0.20 


90 Zr 


fi/A (MeV) 


8.71 


8.69 


8.68 




fich (fm) 


4.26 


4.26 


4.25 




Rn-Rp (fm) 


— 


0.11 


0.09 


116 Sn 


fi/A (MeV) 


8.52 


8.48 


8.50 




fich (fm) 


4.63 


4.60 


4.60 




Rn-Rp (fm) 


— 


0.17 


0.13 


132 Sn 


(MeV) 


8.36 


8.37 


8.34 




-Rch (fm) 




4.70 


4.71 




Rn-Rp (fm) 




0.35 


0.27 


208p b 


(MeV) 


7.87 


7.88 


7.89 




R ch (fm) 


5.50 


5.51 


5.52 




Rn-Rp (fm) 




0.28 


0.21 



TABLE II: Experimental data for the binding energy per nucleon and the charge radii for the 
magic nuclei used in the least square fitting procedure. In addition, predictions are displayed for 
the neutron skin of these nuclei. 



rather by their linear response. Indeed, as argued in Refs. 23, 23( the success of the NL3 set 
in reproducing the breathing mode in 208 Pb is accidental, as it results from a combination 
of both a stiff equation of state for symmetric nuclear matter and a stiff symmetry energy. 
That this is indeed the case may be appreciated by examining the breathing mode in 90 Zr 
- a nucleus with well-developed monopole strength but rather insensitive to the symmetry 
energy — and the energy of the isovector giant dipole resonance (IVGDR) in 208 Pb — an 
observable sensitive to the density dependence of the symmetry energy. In Table 11111 rel- 
ativistic random phase approximation (RPA) results for the centroid energies of the giant 
monopole resonance (GMR) in 208 Pb and 90 Zr, as well as the IVGDR (peak energy) in 208 Pb 
are reported. These small-amplitude modes represent the linear response of the mean field 
ground state to a variety of probes [i^, 24]|. As alluded earlier, the NL3 parameter set 



predicting both a stiff equation of state for symmetric nuclear matter (K = 11\ MeV) and a 
stiff symmetry energy R n — R p = 0.28 fm — overpredicts the centroid energy of the GMR in 
90 Zr and underpredicts the peak position of the IVGDR in in 208 Pb. (Note that a symmetry 
energy that rises rapidly with energy predicts low values for the symmetry energy at the 
densities probed by the isovector dipole mode). In contrast, the good agreement between 
FSUGold and experiment is due to the addition of the two extra parameters (£ to reduce the 
value of K and A v to soften the symmetry energy). Note that an additional softening of the 
symmetry energy may further improve the agreement with experiment. Thus, our present 
prediction of R n — R p = 0.21 fm could be regarded as an upper bound. This smaller value 
for the neutron skin in 208 Pb, generated from the softer symmetry energy, is significant as it 
brings covariant meson-baryon models closer to nonrelativistic predictions based on Skyrme 
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FIG. 1: Proton (charge) and neutron (point) densities for 208 Pb using a variety of values for the 
isoscalar-isovector coupling constant A v . 



Nucleus 


Observable 


Experiment 


NL3 


FSUGold 


208pk 
90 Zr 

208 Pb 


GMR (MeV) 
GMR (MeV) 
IVGDR (MeV) 


14.17 ±0.28 
17.89 ± 0.20 
13.30 ±0.10 


14.32 
18.62 
12.70 


14.04 
17.98 
13.07 



TABLE III: Centroid energies for the breathing mode in 208 Pb and 90 Zr, and the peak energy for 
the IVGDR in 208 Pb. Experimental data are extracted from Refs. 

and 0. 



parameterizations [2l|. The Parity Radius Experiment (PREX) at the Jefferson Laboratory 
should provide a unique observational constraint on the density dependence of the symmetry 
energy. 

Having constructed the accurately calibrated FSUGold parameter set with an equation 
of state (EoS) that is considerably softer than NL3, we now examine some of its predictions 
for a few neutron-star observables. The structure of spherically-symmetric neutron stars 
in hydrostatic equilibrium is described by a solution to the Tolman-Oppenheimer-Volkoff 
(TOV) equation, that may be expressed as a coupled set of first-order differential equations 
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of the following form: 



— = -G 



dP 

dr 
dM 



£(r)M(r) 



P(r) 
S(r) 



1 + 



47ir 3 P(r) 
M(r) 



2GM(r) 



dr 



Anr 2 £(r) . 



(2a) 
(2b) 



In the above equations G is Newton's universal gravitational constant, while P(r), £{r), 
and M(r) represent the pressure, energy density, and enclosed-mass profiles of the star, 
respectively. The second of the two equations [Eq. ()2b[) ] is the differential form for the 
definition of the mass enclosed M(r) up to a radius r. That is, 



dM 
dr 



47rr 2 £(r) 



M(r) 



4:7rx 2 £(x)dx . 



(3) 



The dynamical information is contained in the second of the two equations [Eq. (|2a|)] that 
embodies the notion of hydrostatic equilibrium, namely, the gravitational attraction on any 
mass element of the star must be compensated exactly by the pressure gradient (dP/dr <0). 
The last three terms on the right-hand side of this equation (enclosed in brackets) incorporate 
corrections from general relativity on to the simple Newtonian dynamics [25] . Incorporating 
these three correction terms is critical, as typical escape velocities from neutron stars are of 
the order of half the speed of light. That is, 



^csc/c 



3M 
~R~ 



1/2 



(4) 



where in the above expression the mass of the star (M) is measured in solar masses and the 
radius (R) in kilometers. 

The TOV equations, together with their associated boundary conditions, i.e., P(r = 0) = 
P c and M(r = 0) = 0, are incomplete without an equation of state P = P{£) to specify the 
relation between the pressure and the energy density of the system. Thus, various properties 
of neutron stars, such as their masses, radii, and temperature, probe directly the dynamics 
of neutron-rich matter over a wide range of density. Thus, if general relativity is assumed 
valid — a very modest and safe assumption — then the structure of a neutron star is solely 
determined by the EoS of neutron-rich matter in beta equilibrium. 

For the uniform liquid phase in the mantle (or outer core) we assume an EoS consist- 
ing of neutrons, protons, electrons, and muons in beta equilibrium. Moreover, we assume 
that this description remains valid in the high-density inner core. Thus, transitions to 
exotic phases, such as meson condensates, hyperonic matter, and/or quark matter, are 
are not considered here. Yet at the lower densities of the inner crust the uniform sys- 
tem becomes unstable against density fluctuations, as it becomes energetically favorable 
to separate the system into regions of high- and low-density matter. In this nonuniform 
region the system may consists of a variety of complex structures, such as spherical, cylin- 



drical, rods, plates, etc — collectively known as nuclear vasta 1261 127L While microscopic 



calculations of the nuclear pasta are now becoming available |28l . I29L I30L |31|, it is prema- 
ture to incorporate them in our calculation. Hence, following the procedure adopted in 
Ref. 13211 - a simple polytropic equation of state is used to interpolate between the outer 
crust [33j and the uniform liquid. For an accurate rendition of the expected structure 
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FIG. 2: State-of-the-art rendition of the structure of a neutron star. 

of a neutron star see Fig. |21 (This figure, courtesy of Dany Page, may be found at 
http : //www. astroscu.unam.mx/neutrones/NS-Picture/NS-Picture .html ). 

In Fig. El we display the equation of state for symmetric nuclear matter (left panel) and 
the symmetry energy (right panel) for the NL3 and FSUGold parameter sets. Note that the 
equation of state for pure neutron matter is to an excellent approximation equal to the sum 
of the two: 

^pnm(p) ^ E SNM (p) + S(p) . (5) 

Symmetric nuclear matter saturates for both parameter sets at a baryon density of po — 
0.148 fm -3 and a binding energy per nucleon of E/A— M = — 16.2 MeV. Yet the addition of a 
quartic vector-meson coupling ( yields a significant softening of the equation of state. Indeed, 
the nuclear matter incompressibility gets reduced from the NL3 value of K — 271 MeV to 
i^ = 230 MeV for the FSUGold set. As we shall see, differences in the predictions of various 
neutron-star properties between these two models are also significant. 

Large differences in the density dependence of the symmetry energy are also predicted by 
these two models. An often used parametrization of the density dependence of the symmetry 
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Symmetric Nuclear Matter Symmetry Energy 




FIG. 3: Equation of state for symmetric nuclear matter (left panel) and the density dependence 
of the symmetry energy (right panel) for the NL3 and FSUGold parameter sets. The inset on the 
right-hand panel is a fit over the low-density region of the predicted density dependence of the two 
models; see text for details. 



energy is given by 3J, 35|, 2s 37 , |38[ : 



, 7 
P " 



where Sq is the value of the symmetry energy per nucleon at saturation density and gamma 
is the parameter that quantifies its density dependence. A fit to the density dependence 
of the symmetry energy over the < p < po region yields the following values for the two 
models: 

f36.9MeV forNL3, , fo.98 forNL3, 

on = < and 7 = < 7 

[32.7 MeV for FSUGold, ' [0.64 for FSUGold. V ; 

As alluded earlier, the larger the value of 7 the larger the neutron skin of 208 Pb. Indeed, an 
empirical fit to a large number of mean- field models yields (39^ : 

R n - R p * (0.227 + 0.06) fm * <{' - 28 fm f ° r NL3 ' (8) 

' 0.21 fm for FSUGold. V ; 

At low densities the uniform system becomes unstable against density fluctuations. That 
is, at low densities it becomes energetically favorable for the uniform system to separate 



8 



into regions of high and low densities. Results for the transition density from the uniform 
mantle to the nonuniform crust as a function of the neutron skin in 208 Pb are displayed 
in Fig |3J Various models are used to show the nearly model-independent relation between 
these two seemingly distinct observables. The figure displays an inverse correlation between 
the neutron-skin and the transition density found in Ref 0]. This correlation suggests 
that models with a stiff equation of state (such as NL3) predict a low transition density, 
as it becomes energetically unfavorable to separate nuclear matter into regions of high and 
low densities. Finally, this "data-to-data" relation illustrates how an accurate and model- 
independent determination of the neutron skin of 208 Pb at the Jefferson Laboratory - 
assumed here purely on theoretical biases to be R n — R p = 0.20 fm — would determine an 
important neutron-star observable. 




R -R (fm) 

n p v ' 

FIG. 4: Transition density from the uniform mantle to the non-uniform crust as a function of the 
neutron skin of 208 Pb. Different models are used to show the largely model independent relation 
between these two quantities. 

A mechanism that may also benefit from an accurate determination of the neutron radius 
in 208 Pb is the cooling of a neutron star. Proto-neutron stars are created hot in supernova 
explosions but then cool rapidly through neutrino emission |4ol |. Indeed, 99% of the energy 
released in a supernovae explosion is carried away by neutrinos. That proto-neutron stars 
are born with very high temperatures was inferred from the few detected neutrinos from 
SN1987A that suggest a neutrinosphere temperature as high as 5 MeV ^jj. Initially, the 
proto-neutron stars will cool rapidly via the direct URCA process — an efficient cooling 
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mechanism that consists of neutron beta decay followed by electron capture 42l |: 



n — > p + e + v e , (9a) 
p + e~ — > n + u e . (9b) 

As the proto-neutron star material becomes neutron rich — a process known as "neutron- 
ization" — the direct URCA process ceases and is replaced by the modified URCA reaction: 

n + n — > n + p + e~ + v e . (10) 

The modified URCA process, however, is relatively slow as a second nucleon is necessary 
to conserve both energy and momentum at the Fermi surface |43 |. Incidentally, the term 
URCA was coined by an Ukrainian — George Gamow. After Mario Schoenberg and George 
Gamow went gambling at the now defunct URCA Casino in Rio de Janeiro, Gamow was 
so impressed by the roulette table where money just disappeared that he said: "well, the 
energy disappears in the nucleus of the supernova as quickly as the money disappears at 
that roulette table". I have also been told that in Russian slang URCA can also mean a 
pickpocket — an individual that can steal your money in a matter of seconds. 

Recent X-ray observations of the neutron star in 3C58 [44|, Vela j45[, and Geminga [I| 



indicate low surface temperatures. Moreover, the low quiescent luminosity in the transiently 
accreting binaries KS 1731-260 and Cen X-4 |48[ suggest rapid cooling. As X-ray 
observatories progress and our knowledge of neutron-star atmospheres and ages improve, 
additional "cold" neutron stars may be discovered. Such low surface temperatures appear 
to require enhanced cooling from reactions that proceed faster than the modified URCA 
process of Eq. (fTU|) . 

Enhanced cooling may occur via the weak decay of additional hadrons such as pion or 
kaon condensates, hyperons, and/or quark matter [i]J 13] . Yet the most conservative 
enhanced-cooling scenario is the direct URCA process of Eq. 0. This mechanism is not 
"exotic" as it only requires protons, neutrons, electrons, and muons — constituents known to 
be present in dense matter. However, to conserve energy and momentum at the Fermi surface 
the sum of the Fermi momenta of the protons plus that of the electrons must be greater than 
the neutron Fermi momentum. This requires a relatively large proton fraction Y p = Z/A. 
As the symmetry energy represents a penalty imposed on the system for departing from the 
symmetric N = Z limit, the proton fraction is highly sensitive to the density dependence 
of the symmetry energy. In essence, the stiff er the symmetry energy the larger the proton 
fraction — as it is energetically unfavorable for the proton fraction to depart significantly 
from its symmetric Y p = l/2 value. 

In Fig. El we display the proton fraction Y p as a function of baryon density for the NL3 
and FSUGold parameter sets. The location of the symbols ("stars") and the quantities 
enclosed in parenthesis indicate the central density and the minimum value of the proton 
fraction necessary for the direct URCA process to be effective in a M = 1.4M Q neutron 
star. In particular, if the solid curve passes above the symbol, then the direct URCA 
process is allowed in such a neutron star. Several features of this plot are worth focusing 
on. First, because of its considerable stiffer symmetry energy, the proton fraction predicted 
by the NL3 set increases much more rapidly than for the FSUGold model. Second, not only 
does the symmetry energy increases more rapidly in the case of the NL3 set, but so does 
the equation of state for symmetric nuclear matter (see Fig. EJ). As a result, the pressure 
required to support a M = 1.4M neutron star against gravitational collapse is reached at 
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FIG. 5: Proton fraction as a function of baryon density for the NL3 and FSUGold parameter sets. 
The "stars" indicate the central density and the minimum value of the proton fraction necessary 
for the direct URCA process to occur in a 1.4 solar- mass neutron star. 



a much higher density in the FSUGold model than in the NL3 set: p c = 3.6 p f° r FSUGold 
vs p c = 1.9 po f° r NL3. It is precisely this fact — the much higher density reached in the 
case of the FSUGold parameter set — that allows a M = 1.4M to be cooled by the direct 
URCA process in spite of its soft symmetry energy. Yet, either a softer equation of state for 
symmetric nuclear matter or a softer symmetry energy may preclude the enhanced cooling 
of a M = 1.4M Q neutron star by the (nucleon) direct URCA process. If such is the case, then 
the enhanced cooling of a M = 1.4M Q neutron star may provide strong evidence in favor of 
exotic matter (e.g., "quark matter") in the core of a neutron star. We close this section by 
listing in Table ITVl the predictions of both models for a variety of neutron-star observables. 



IV. CONCLUSIONS 



In conclusion, a new accurately calibrated relativistic model ( "FSUGold" ) has been fitted 
to the binding energies and charge radii of a variety of magic nuclei. In this regard, the new 
parametrization is as successful as the NL3 set which has been used here as a useful paradigm. 
In particular, symmetric nuclear matter saturates at a Fermi momentum of fcp = l-30 fm _1 
(corresponding to a baryon density of 0.15 fm~ 3 ) with a binding energy per nucleon of 
B/A=— 16.30 MeV. Further, by constraining the FSUGold parameter set by a few nuclear 
collective modes, we obtain a nuclear-matter incompressibility of K = 230 MeV and a neutron 
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Neutron-Star Observable 


NL3 


FSUGold 


p c (fm- 3 ) 


0.052 


0.076 


R (km) 


15.05 


12.66 


M max (M ) 


2.78 


1.72 


PuRCA ( fm_3 ) 


0.21 


0.47 


M URC a(M ) 


0.84 


1.30 


AMtjrca 


0.38 


0.06 



TABLE IV: Predictions for a few neutron-star observables. The various quantities are as follows: 
p c is the transition density from nonuniform to uniform neutron-rich matter matter, R is the radius 
of a 1.4 solar-mass neutron star, M max is the limiting mass, p URCA is the threshold density for the 
direct URCA process, Murca is the minimum mass neutron star that may cool down by the direct 
URCA process, and AMtjrca is the mass fraction of a 1.4 solar-mass neutron star that supports 
enhanced cooling by the direct URCA process. 

skin thickness in 208 Pb of R n — R p = 0.21 fm. While the description of the various collective 
modes imposes additional constraints on the EoS at densities around saturation density, 
the high-density component of the EoS remains largely unconstrained. Thus, we made no 
attempts at constraining the EoS at the supranuclear densities of relevance to neutron-star 
physics. Rather, we simply explored the consequences of the new parametrization on a 
variety of neutron star observables and eagerly await high-quality data that will constrain 
the high-density component of the EoS. In particular, we found a limiting neutron-star mass 
of M max = 1.72M Q , a radius of R— 12.66 km for a M = 1.4M Q neutron star, and no direct 
URCA cooling in neutron stars with masses below M = 1.3M Q . It is interesting to note 
that recent observations of pulsar-white dwarf binaries at the Arecibo observatory suggest 
a pulsar mass for PSRJ0751+1807 of M = 2.1+^M at a 95% confidence level Q. If this 
observation could be refined, not only would it redefine the high-density behavior of this 
(and many other) EoS, but it could provide us with a precious boost in our quest for the 
equation of state. 
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